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I. INTRODUCTION 

Many proposals for explaining the apparent shortcom- 
ings of the Standard Model have been advanced. The 
Unparticle Model proposed by Georgi [l[ aimed to in- 
clude in the Standard Model massive but scale-invariant 
particles, sharing the same physics of the scale-dependent 
counterparts. These objects, called 'unparticles', could 
play an important role in low-energy physics Q , since the 
model implies that unparticles can be exchanged between 
massive particles, leading to a new force called 'ungrav- 
ity'. This "fifth" force would add a perturbation term to 
the newtonian gravitational potential, although the ex- 
act potential can not be obtained because the distance 
at which the perturbed potential matches the newtonian 
expression needs to be known. In order to bypass this 
limitation, the perturbed potential has been assumed to 
be of the form H 
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where djj (the scaling dimension of the unparticles op- 
erator Ojj) is ~ 1, as a reasonable approximation, and 
Rg is the characteristic length scale of ungravity, given 
by 



A/* with du — 1- For this purpose, a suitable quasi- 
newtonian gravitational system needs to be studied and 
compared with the pure newtonian results. A first study 
of this regime by Bertolami, Paramos and Santos 3] has 
addressed the stellar equilibrium problem, deriving a per- 
turbed Lane-Emden equation further applied to the Sun. 
They explored the well-known similarity of the full stellar 
structure to an n = 3 polytropic model, and derived lim- 
its from the maximum allowed uncertainty in the central 
temperature AT C /T C = 0.06. 

In spite of the successful derivation of meaningful lim- 
its to the unparticle parameters, it is known that the 
detailed structure of the Sun is actually quite compli- 
cated, and many physical factors have to be considered 
beyond the simplest Chandrasekhar's polytropic model 
[j]. Therefore, it is worth considering another very well- 
known system to which the Chandrasekhar theory gives 
an even better representation: the white dwarf sequence. 
We shall show below that an important feature of these 
sequences (the maximum mass) is sensitive to the un- 
particle quantities and allows to impose strong limits on 
them. 



II. STELLAR EQUILIBRIUM AND 
WHITE-DWARFS 
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where Ajj is the energy scale of the unparticle interac- 
tions, Mpi = 2.4 x 10 18 GeV is the Planck mass and a is 
a constant dependent on the type of propagator consid- 
ered. 

The problem addressed in this work is to determine 
the bounds of the mass of the interaction (un)particle 



Since Chandrasekhar's polytropic model is widely 
known, we briefly recall how the unparticle theory mod- 
ifies it, as first shown by Bertolami, Paramos and Santos 
Q. The equations of stellar hydrostatic equilibrium and 
mass conservation can be reduced to a second order dif- 
ferential equation, if a polytropic equation of state of the 
form P = Kp 1+1 / n is assumed to hold. If the density is 
written as p = p c n , and the radius as r — /3£, one can 
easily obtain the original Lane-Emden equation 
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where n is the polytropic index, £ is the dimensionless 
radius and f3 is given by 
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K being the polytropic constant dependent on the spe- 
cific value of n. This equation is subject to the usual 
boundary conditions: p(r = 0) = p c and dP/dr = 
for r — 0, which translates to 0(£ = 0) = 1 and 
0'(£ = 0) = d9/d£_ = 0. The detailed derivation of the 
Lane-Emden equation can be consulted in the classical 
reference [|J. These results can be used to derive the 
mass-radius relation, given by Eq.©. 
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Considering that white-dwarfs are small stars com- 
posed by electron-degenerate matter, in which the core 
material no longer undergoes fusion reactions, the Lane- 
Emden equation describes very well their behavior as a 
result of the proximity of the electronic component to a 
polytropic form. As is well-known, in the non-relativistic 
limit the white dwarf matter can be represented by a 
n = 3/2 polytrope and in the relativistic case the n = 3 
is quite accurate. 

The same technique employed by Chandrasekhar can 
be used to obtain the perturbed Lane-Emden equation. 
From the perturbed potential given by Eq.©, the per- 
turbed gravitational acceleration can be obtained via 
F = — VV @, used to determine the hydrostatic equilib- 
rium equation, as seen in Eq.©. 
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Then, from this hydrostatic equilibrium equation, 
Bertolami, Paramos and Santos derived a perturbed 
Lane-Emden equation 
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in order to examine meaningful bounds on the parame- 
ters Rg in Eq. ([2]) , from which one can obtain the bounds 
of the mass of the interaction (un)particle M* based on 
astrophysical constraints (the parameter Zg = Rg /ft has 
been defined here to simplify the expressions). We have 
corrected a slight misprint in the paper Q which would 
preclude the derivation of Eq.©, related to the definition 
of the variable /3 as given by © . 



Once these matters are settled, the perturbed mass- 
radius relation can be obtained in the usual manner, be- 
ing formally identical to the analogous mass-radius re- 
lation given by Chandrasekhar (Eq.©), but with the 
perturbative effects present through the solutions of the 
perturbed Lane-Emden equation and its derivative 
and 0'(£i)y, respectively, both evaluated at the first zero 
as in Ref. [J] . This modification sets the stage for an anal- 
ysis leading to novel bounds on the unparticle parame- 
ters. Since we are interested on the maximum possible 
mass for a white dwarf, we set the parameters related 
to the Lane-Emden model to the relativistic limit (poly- 
tropic index n = 3), and the white dwarf mass-radius 
relation becomes 



(8) 



losing any dependence on the radius, as expected. 

The key feature pointed out above, related to the ex- 
plicit dependence of the maximum mass through the zero 
of the Lane-Emden function and its slope, leads to quite 
a strong dependence of the maximum mass with djj and 
Rg, which can be substantially different from the "canon- 
ical" M Ch = 1.457(2/ Me ) 2 M 0,0. 



III. OBSERVATIONS AND ANALYSIS 

Following the framework presented by Bertolami, 
Paramos and Santos , we obtained the numerical solu- 
tions into Eq.© varying both parameters Zg and du, and 
used these solutions on Eq.© to determine the masses 
associated to the combination of du and Zg- It is impor- 
tant to note that the combination Zg and djj generated 
a broad range of masses, spanning from M « 1.291 M Q 
to M « 1.874 Mq. This in turn means that it is possi- 
ble to constrain the parameters based on the maximum 
mass for white dwarfs, because this mass must complain 
with the values obtained by observational data. Follow- 
ing this reasoning, we assumed three reasonable values 
for the maximum mass 

• the maximum observed masses in a large white 
dwarf sample, from the recent work by Kepler and 
collaborators 0], is M = 1.33 M . Neverthe- 
less, in order to effectively constrain the values of 
djj near 1, we had to consider masses of at least 
M = 1.36 Mq, otherwise djj would be far from 1 
and the gravitational corrections in Eq.© would 
be too large; 

• the "canonical" limit, Men — 1-457 Mq for carbon- 
type white dwarfs; in spite of this widely accepted 
value is still beyond the actual observed maximum, 
its use seems very reasonable 
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• an even larger value, above the former "Chan- 
drasekhar's" limit, M — 1.60 Mq, arbitrarily cho- 
sen to represent an extreme limit allowed in nature. 
If true, the number of objects between the actually 
observed maximum M — 1.33 Mq and this pro- 
posed extreme value of M = 1.60 M© must be 
substantial in a large sample such as the one ana- 
lyzed by Kepler et al. [7j, although none has been 
actually reported. 

From the data above mentioned, we constructed the 
contour plots where the combination of dy and £g gen- 
erates the desired maximum mass. This plots are de- 
picted in FigQ] for ranges larger than the radius of the 
star (£g > ^, djj < 1) and in Fig [3] for ranges smaller 
the radius of the star (£q < 1, du > 1). One should 
note that for djj > 1, only masses of M = 1.36 Mq 
and M = 1.457 M© are obtained, while a mass of 
M = 1.6 Mq can be imposed for du < 1 only. 
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FIG. 1: Contour plot (locus) of maximum mass (M = 
1.6 Mq) in the Rc-du plane, for an interaction range larger 
than the radius of the star. The dashed area gives masses 
bigger than M = 1.6 Mq, and should be considered as for- 
bidden. 

The connection with the (un)particle mass follows by 
solving equation Eq.©, for the ratio M^/Mpi, namely 
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which, by using Rc(du) defined in the plots above, it is 
possible to infer the mass of the interaction (un)particle 
M* . We point out that this procedure differs significantly 
from the method used by Bertolami, Paramos and Santos 
Q, where the 6% uncertainty leads to terms R- (du) and 
R+(du), because we use exact values for the maximum 
masses. Nevertheless, it is clear that the ratio must be 
interpreted as a lower bound to the ratio M^/Mpi, con- 
sidering that lower maximum masses (which would be in 
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FIG. 2: Contour plot (locus) of maximum masses M = 
1.36 M@ and M = 1.457 Mq in the RG-du, for an inter- 
action range smaller than the radius of the star. The curve 
spanning from 1.055 < du < 1.06 (upper right corner of the 
graphic) produces a maximum mass of M = 1.36 Mq and 
the curve between 1 < du < 1.06 (lower side of the graphic) 
the mass of M = 1.457 Mq. The upper white area gives 
maximum masses lower than M = 1.36 Mq, the dashed area 
yields masses in the range 1.36 Mq < M < 1.457 Mq, and 
the lower area produces masses bigger than 1.457 Mq. If one 
accepts the latter value as the maximum allowed, the dark 
area is forbidden 



conflict with observations, as stated before) would result 
in even lower ratios. Plotting the ratio as a function of 
du and fixing the parameters a and Au in the same way 
as Bertolami, Paramos and Santos, we obtain the ranges 
for the mass of the interaction (un)particle depicted in 
FigEJand FigfJ]for the case of the "canonical" mass and 
two different choices of Au- 




FIG. 3: Lower bounds on log(Af,/Af p i), for a = 0, a = 2/3, 
a = 1, a = 1.9, from top to bottom respectively. The fixed 



values are Au = 1 TeV and M = 1.457 M, 
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Analyzing the data for M = 1.36 Mq, dus fall in the 
range 1.055 < du < 1.060. Therefore, this maximum 
mass yields a lower bound to M* in the range (0.03 — 
0.50)M p; . For M = 1.457 Mq, d v s lies in the range 
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FIG. 4: Lower bounds on log(M»/M p! ), for a = 0, a = 2/3, 
a = 1, a = 1.9 from top to bottom respectively. The fixed 
values are Ku = 10 3 TeV and M = 1.457 M Q . 

1 < djj < 1.06, yielding a lower bound to M* in the range 
(0.1 — 1.6)M p i. The last case, with an assumed maximum 
mass M = 1.6 M© and the interaction range larger than 
the star radius, the only allowed values for djjs are in the 
range 0.94 < djj < 0.982, and lower bounds to M* in the 
range (0.5 - 6.3)M pi . 

IV. CONCLUSIONS 

We have shown in this work that quite strong limits 
to the unparticle parameters can be obtained by using a 
simple form of the polytropic theory of Chandrasekhar 
adding a perturbation to the Lane-Emden equation, as 
first obtained by Bertolami, Paramos and Santos Q, and 
applying it to the white dwarf sequences. 



The key point elaborated here is that a change on the 
unparticle parameters would affect the maximum mass 
allowed to white dwarfs, and thereby we explored this 
characteristic in order to limit the values of such param- 
eters. 

The requirement that the maximum mass can not be 
too small (because it would conflict with a few massive 
stars [7[ or too big (because it would lead to unobserved 
supermassive white dwarfs) limit the values of M* to a 
confidence range of O.lMp; < M* < l.QMpi from this 
analysis alone for the case djj > 1. For the case dxj < 1, 
the mass-radius relation gives only masses bigger than 
the canonical value. Considering that until today there 
is no observation of white dwarfs with such high masses, 
this analysis may be interpreted to mean that values of 
d\j < 1 are not allowed. 

Following a different approach, based on a cosmological 
scenario, Bertolami and Santos [l3[ considered the varia- 
tion of the gravitational coupling at the time of big bang 
nucleosynthesis, tensor exchange and the scaling dimen- 
sion djj = 1.1 and found M* to be > 0.05Afp;, which is 
very close to the bounds found for 1 < dy < 1.06. Other 
works studying complementary bounds [!, [||-[l2[ could be 
combined to address the viability of a general unparticle 
model, unless one can manage to evade the bounds al- 
together. Even if so, a general argument to constrain 
the admissible perturbations to the newtonian potential 
can be made via the perturbed Lane-Emden equation, 
resorting to the observed massive white dwarfs. 
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